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Four-point measurements using a multi-tip scanning tunneling microscope (STM) are carried out 
in order to determine surface and step conductivities on Si(lll) surfaces. In a first step, distance- 
dependent four-point measurements in the linear configuration are used in combination with an 
analytical three-layer model for charge transport to disentangle the 2D surface conductivity from 
non-surface contributions. A termination of the Si(lll) surface with either Bi or H results in the 
two limiting cases of a pure 2D or 3D conductance, respectively. In order to further disentangle the 
surface conductivity of the step-free surface from the contribution due to atomic steps, a square four- 
probe configuration is applied as function of the rotation angle. In total this combined approach leads 
to an atomic step conductivity of (Jstep = (29 ± 9) and to a step-free surface conductivity 

of Usurf = (9 ± 2) • 10“® for the Si(lll)-(7x7) surface. 


The increasing importance of surface conductance 
compared to conductance through the bulk in modern 
nanoelectronic devices calls for a reliable determination 
of the surface conductivity in order to minimize the in¬ 
fluence of undesired leakage currents on the device per¬ 
formance or to use surfaces as functional units. A model 
system for corresponding investigations is the Si(lll)- 
(7x7) surface. Over the years a wide range of values 
for the conductivity of this surface has been reported, 
spanning several orders of magnitude [T], and the latest 
measurements deviate still by a factor of 2 to 3 nig. 
The main difficulty in measuring the surface conductiv¬ 
ity is to separate the 2D conductance at the surface from 
the conductance through other channels, e.g. the bulk 
and the space charge layer. 

Here, we use a four-tip scanning tunneling microscope 
[1] for distance-dependent measurements of the four- 
point resistance on Si(lll), as shown in the inset in Fig.[^ 
for a linear tip arrangement, in combination with a three- 
layer model for charge transport. This method allows the 
separation of the surface conductance from other contri¬ 
butions due to the characteristic probe spacing depen¬ 
dency of different conductance channels. Further on, we 
analyze the anisotropy of the surface conductance caused 
by the influence of atomic steps, which allows to deter¬ 
mine the conductivity of a single step and the step-free 
surface. 

Analytic equations relating the measured four-point re¬ 
sistance to a conductivity can be obtained easily for pure 
2D or 3D geometries, i.e. four tips positioned on a con¬ 
ducting sheet (surface) or on a half space (bulk), as [1] 

and R% = ^ -s-i (1) 

7r(T2n 2Tra3D 

with an equidistant probe spacing s, the 2D surface con¬ 
ductivity <J 2 d and the 3D bulk conductivity o^d- The 
equation for the 2D case shows the hallmark of a 2D 


channel, namely the fact that the surface conductance is 
independent of the probe spacing, while the conductance 
through a 3D channel depends on the distance s of the 
four probes. In order to minimize the number of tips to 
be repositioned, we preferentially use a non-equidistant 
spacing, in which three tips remain at a mutual distance 
of s = 50 /xm, while only the distance x between tip 1 and 
tip 2 is varied (Fig. [^. In this non-equidistant setup the 
hallmark of the constant four-probe resistance is lost for 
the 2D case, since Eq. [^has to be modified as shown in 
[5H7] and summarized in the Supplemental Material [5]. 

The four-point resistance measured on a Si(lll)- 
(•\/3 X •\/3)R.30° Bi-terminated (1 ML) surface of an 



FIG. 1. (Color online) Four-point resistance of a Bi- 
terminated Si(lll)-(%/3 X -/3)R30° sample as fnnction of the 
probe distances s and x for the eqnidistant (right half) and 
the non-equidistant configuration (left half). The red solid 
line represents the behavior expected for a pure 2D conduc¬ 
tance with (TBi = (1.4 ± 0.1) • 10“^D“^/n. In the inset the 
linear measurement configuration is shown. 
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n-doped sample (2kl7cm) is shown in Fig. for the 
non-equidistant configuration with distances x < s = 
50 fim, and for the equidistant configuration with dis¬ 
tances x = s > 50 /rm (details of sample preparation 
and measurement procedure are described in the Sup¬ 
plemental Material 0 ). The constant behavior in the 
equidistant range s > 50 /rm indicates a pure 2D char¬ 
acter of conductance. Annother indicator for 2D sur¬ 
face transport is the fact that the four-point resistance, 
which is expected considering only the bulk conductivity, 
is several orders of magnitude larger than the observed 
one. Therefore, we compare the experimental data to 
a 2D model, and a good correspondence is obtained for 
CBi = (1.4±0.1) • lO^'^D^^/n (solid red line) confirming 
that the charge transport in the Bi-terminated Si(lll) 
sample occurs almost exclusively through the 2D surface 
channel. Similar results were found for two differently 
doped samples. 

Subsequently, the distance dependence of the four- 
point resistance was measured on a clean Si(lll)-(7x7) 
sample. The results for an n-doped sample (700ncm) 
are shown in Fig. The observed decreasing four-point 
resistance for increasing equidistant probe spacing s indi¬ 
cates that a non-surface channel contributes to the charge 
transport, since a pure 2D conduction exhibits a constant 
behavior in the equidistant region (cf. Fig.[^. Thus, the 
measured four-point resistance should be modeled by a 
conductance channel through the surface states as well as 
additional contributions from the bulk and a near-surface 
space charge layer. However, in this case Eq. cannot 
be applied. 

Often an approximation of a parallel circuit consisting 
of the four-point resistance of the surface and the bulk 
(plus space charge layer) is used [B], but this approach 
has two shortcomings. First, a complete separation of the 
surface conduction channel and the bulk is assumed. Sec¬ 
ond, the two-point resistances, not the four-point resis¬ 
tances, determine, which amount of current flows through 
the surface layer and which part through the bulk/space 
charge layer. So, the preferred way for the current to split 
up depends on the details of the injection, e.g. the size 
of the current injecting contact (tip diameter) [3]. Thus, 
if more than one current path exists, the four-point re¬ 
sistance depends on possible transitions between charge 
transport channels as well as on the properties of the cur¬ 
rent injecting contacts, so that the well-known statement 
that the four-point resistance measured on the surface is 
independent of the contact resistances is not completely 
true. 

In order to describe the charge transport through the 
different channels more accurately, we use an analytically 
derived three-layer model for conductance. In this model, 
the bulk enters with its known conductivity, while the 
surface conductivity is the parameter to be determined 
by a ht to the data. The space charge layer is approx¬ 
imated by an intermediate layer with a certain thick- 



FIG. 2. (Color online) Four-point resistance of an n-doped 
Si(lll)-(7x7) sample as fnnction of the probe distances s and 
X for the eqnidistant and the non-eqnidistant confignration. 
A three-layer model for charge transport yields the solid red 
line with (T 2 d = (5.1 ± 0.7) • 10“® located between the 

two limiting cases of a pnre surface conductance (dashed bine 
line) and a pure bulk conductance (dashed green line). If the 
Si(lll) surface is hydrogen-terminated, the surface contribu¬ 
tion vanishes completely and the remaining bulk conductance 
can directly be measured (gray data points). 

ness and conductivity. These two values are obtained by 
the solution of the Poisson equation, which considers the 
known Fermi level pinning of the Si(lll)-(7x7) surface 
mm- Nevertheless, the use of the bulk doping concen¬ 
tration as an initial parameter in this calculation in order 
to determine the width and average conductivity of the 
space charge layer is not sufficient to describe the data in 
Fig.H However, it is known that high-temperature an¬ 
nealing up to 1200 °C performed for cleaning the Si(lll) 
surface causes a dopant redistribution and an additional 
p-type doping in the near-surface region due to boron 
in-diffusion [HHE] or possible formation of near-surface 
single vacancies m- These effects can lead to a reduced 
carrier concentration in the space charge layer. Gener¬ 
ally, the details of the modification of the near-surface 
doping depend highly on the specific method and setup 
used for sample preparation. We find that the experi¬ 
mental data can be described well for a conductivity of 
2.5 • 10“^ and a thickness of 3.1 /rm for the inter¬ 

mediate layer representing the space charge layer. This 
quite approximate modeling of the space charge region 
as only one layer with constant conductivity seems to be 
sufficient, as the surface conductivity obtained from the 
fit to the measured data turns out to be insensitive to 
the specific properties of the intermediate layer. 

Overall, the three-layer model results in a much more 
accurate description of the measured four-point resis¬ 
tance than the simple parallel-circuit model, since it 
avoids the artificial separation between the surface and 
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the non-surface channels and takes into account the in¬ 
jection geometry giving rise to a charge transport inside 
and between the layers according to their properties. The 
analytical derivation of the model is described in detail 
in the Supplemental Material [S]. 

The best fit to the measured four-point resistance us¬ 
ing the three-layer model is shown as a solid red line in 
Fig-i and results in a surface conductivity of (J 2 D = 
(5.1 ± 0.7) • For comparison the two lim¬ 

iting cases are marked in Fig. The four-point resis¬ 
tance arising from a pure 2D conductivity a^D is shown 
as dashed blue line, while the four-point resistance in¬ 
duced by a pure 3D conductance, with its 1/s behavior 
in the equidistant configuration, is indicated as dashed 
green line featuring a bulk conductivity value, which is 
confirmed by an additional experiment described below. 
In the non-equidistant region the measured four-point re¬ 
sistance is close to the one expected from a pure surface 
conductance (less than 6% deviation for a; < 50/am), but 
for larger probe spacing an increasing deviation from the 
2D behavior is observed. This reflects the well-known 
general tendency that the conductance is more surface- 
dominated for small probe distances, while a non-surface 
contribution develops more significantly for larger dis¬ 
tances [T]. However, the observed four-point resistance 
does not approach the 1/s bulk behavior for s > 50 /im, 
because the space charge layer blocks the charge trans¬ 
port into the bulk due to the low conductivity of the de¬ 
pletion zone. So, the four-point resistance in the equidis¬ 
tant range particularly reflects the properties of the space 
charge layer and the bulk, while the non-equidistant re¬ 
gion is more suitable for the determination of the surface 
conductivity. In total, the three-layer model including 
the intermediate layer describes the experimentally ob¬ 
served behavior very well. Results obtained for other 
doping levels are shown in the Supplemental Material [S] 
and confirm the results presented above. 

An additional experiment is used to explore, if the 
bulk conductivity can be measured directly with the 
four-probe setup after removing the surface conductance 
channel. A hydrogen termination of the Si(lll) surface 
by a treatment in HF, resulting in the formation of the 
Si(lll)-(lx 1)-H, is known to remove the surface states 
present on the 7x7 surface m- The gray data points 
in Fig. show the distance dependence of the four-point 
resistance in the non-equidistant region measured on a 
hydrogen-terminated Si(lll) sample. The dashed green 
line corresponds to a fit using a pure 3D bulk behavior 
with a resistivity of = (580 ± 70) Hem, which is close 
to the macroscopically measured nominal bulk resistiv¬ 
ity of (700±50) Hem and therefore confirms that without 
surface states a pure 3D bulk conductance is obtained. 

While the distance-dependent four-point measure¬ 
ments could disentangle the surface conductivity from 
non-surface contributions to charge transport, the influ¬ 
ence of atomic steps located on the (7x7)-reconstructed 


Si surface has not been considered up to now. The con¬ 
ductivity arising from a single step for a current passing 
through it can be treated as scalar quantity. However, if 
a larger surface area is taken into account, the step array 
leads on average to an anisotropic conductivity described 
by the tensor components crU along the step edges and 
perpendicular to the step edges m- So, the anisotropic 
conductance is a macroscopic (mean field) result of the 
different number of step edges per unit length along dif¬ 
ferent current paths. It turns out that the linear four- 
point measurement configuration (Fig. is not sensi¬ 
tive to a two-dimensional conductance anisotropy m- 
However, in a square arrangement of the four probes, 
as shown in Fig. |^c), an angle-dependent four-point 
resistance is obtained from the solution of the Poisson 
equation for an anisotropic 2D sheet muo] 

-b 1^ - 4cos26»sin26l - l) 

(sin^e -b fF cos 26 »)^ 

with C = l/(47r'\/(Tll(T-*-). 

Results for the measured anisotropic four-point resis¬ 
tance on an n-doped Si(lll)-(7x7) sample (700 Hem) for 
a probe spacing of s = 50 p,m are shown in Fig. [^a) as 
a function of the rotation angle 9 relative to the step di¬ 
rection. The four sets of differently colored data points 
in angle increments of 5° arise from the fact that for 
one fixed orientation of the probes four different rotation 
angles can be realized by successively assigning different 
probes as current and voltage probes. 

A fit of Eq. to the experimental data is shown 
as a dotted blue line in Fig. [^a) describing the an¬ 
gle dependence quite well [H]. Nevertheless, the mean 
field approach applied so far assumes only straight step 
edges. However, the typical step structure present on our 
0.25° misoriented Si(lll)-(7x7) sample surface shown in 
Fig. I^b) consists of steps aligned mainly along two di¬ 
rections with average angles of a « 8° and /3 « 21° with 
respect to the average step orientation (indicated as red 
lines), which now defines trll and a^. This average step 
orientation arises from the macroscopic azimuthal direc¬ 
tion of the sample miscut and is not aligned with the 
low-index orientations of the step edges. To model this 
more complicated non-parallel step pattern, we consider 
as first-order approximation a superposition of two step 
orientations with angles a and /3 relative to the aver¬ 
age step orientation weighted with their respective por¬ 
tion extracted from Fig. ib). This leads to a slightly 
skewed curve shown as solid red line in Fig. [^a), which 
consists of an amount of 70% and 30% of the two sin¬ 
gle contributions, respectively (dotted green lines). The 
model including the two step orientations describes the 
data as well as Eq. but contains a better approxima¬ 
tion of the sample step structure, and results in cr|| = 
(9 ± 2) • 10"® H" Vn and = (1.7 ± 0.4) • 10"® H" VD 


R{9) = C • In 
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FIG. 3. (Color online) (a) Four-point resistance measured on 
a Si(lll)-(7x7) surface in the square configuration (shown in 
(c)) with s = 50/rm as a function of rotation angle 6 between 
the average step orientation and the line connecting the cur¬ 
rent injecting tips (colored data points). The fits to the data 
using either only one average step direction or a superposi¬ 
tion of two step orientations (respective parts are shown as 
green dotted lines) are indicated as dotted blue line and solid 
red line, respectively, (b) STM image of the Si(lll)-(7x7) 
surface showing the representative step arrangement on the 
sample. Two adjacent step edges are highlighted (solid green 
and blue line) consisting of two main step directions indicated 
by the angles a and /3 relative to the average step orientation 
(solid red lines). 


with an anisotropy ratio of cr||/cr_L « 5. The geomet¬ 
ric mean Vcrilcr-*- = (3.9 ± 0.6) • has nearly 

the same value within the error tolerances as the sur¬ 
face conductivity (J 2 d obtained in the linear configura¬ 
tion. Thus, the two independent methods, the distance- 
dependent linear configuration and the angle-dependent 
square configuration, yield the same results for the sur¬ 
face conductivity. 

In a last step, we approximate the measured mean field 
anisotropic conductivity by the scalar resistivities of a 
step-free terrace Psurf and a single step Pstep- Considering 
first the direction parallel to the steps, no step edges have 


to be crossed by the current, which results in the relation 

l/o-ll = pll = Psurf • (3) 

Second, the resistivity perpendicular to the step edges is 
composed of additive contributions from the steps and 
the step-free terraces and can be expressed as series re¬ 
sistance, resulting in [TH] 

l/a^ = p-L = -p p.tep/dj^^p (4) 

with d^ep denoting the average distance between the 
steps. From the two relations in Eq. and Eq. 
finally the conductivity of the step-free Si(lII)-(7x7) 
surface can be disentangled from the influence of the 
step conductivity as (Tgurf = (9 i 2) • and 

CTstep = (29 ± 9) The value of the surface con¬ 

ductivity CTsurf is a factor of 2 to 6 larger than the values 
obtained in recent experiments mm- Such smaller val¬ 
ues may be explained, as these experiments are based 
on a more indirect comparison of the conductivity be¬ 
fore and after quenching the surface states by adsorption 
of atoms/molecules. For the quenched system several 
conditions have to be fulfilled: (a) the surface sates of 
the surface under study are completely quenched, (b) the 
space charge layer conductivity is not influenced by the 
adsorbed layer, and (c) the adsorbed layer induces no 
(additional) surface conductance. If one of these condi¬ 
tions is not fulfilled, these experiments based on the dif¬ 
ference method result in different values for the surface 
conductivity. 

From a comparison of the surface resistivity and the 
step resistivity, the following relation is obtained. The 
resistance of one step (per unit length) corresponds to 
the resistance of a segment of the step-free Si(lll)-(7x7) 
surface (per unit length) of a width of 300 nm. For our 
sample with a step density of 14 steps/pm, the contribu¬ 
tion of the step resistance to the total resistance has a 
substantial amount of 80 % for a current flowing in the 
perpendicular direction. In general, the presence of steps 
will reduce the surface conductivity of the Si(lll)-(7x7) 
considerably, however, in a well predictable manner. 

In conclusion, we combined the distance-dependent lin¬ 
ear configuration for four-point resistance measurements 
on Si(Ill) surfaces with a three-layer model for charge 
transport in order to disentangle the surface conduc¬ 
tivity from non-surface contributions (bulk and space 
charge layer conductivity). The influence of atomic sur¬ 
face steps is obtained by measurements in the angle- 
dependent square configuration resulting in a step-free 
surface resistivity of psurf = (116 ± 26) kfl/D and a step 
resistivity of pstep = (3.4 ± 1) • 10“^ flm for the Si(lll)- 
(7x7) surface. These two generic methods can easily be 
used to determine surface conductivities of other mixed 
2D/3D systems, like different semiconductors or topolog¬ 
ical insulators. 
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In this supplemental material, details about sample preparation and measurement procedure are described, and 
some additional experimental results obtained for differently doped Si samples are presented. Furthermore, the three- 
layer model used to describe the four-point resistance measured in the distance-dependent linear probe configuration 
is discussed. 


DETAILS OF THE SAMPLE PREPARATION AND MEASUREMENT PROCEDURE 


The sample preparations and measurements are carried out under ultra-high vacuum (UHV) conditions with a base 
pressure of ~ 1 • 10“^°mbar. Subsequent to the cleaning process by direct current heating to 1230°C the Si(lll) 
samples are cooled down slowly, especially in the vicinity of the transition temperature of ~ 800 °C, for establishing 
the 7x7 surface reconstruction. Bi passivation is obtained afterwards as described in |Slj to achieve a (-s/S x y3)R30° 
surface terminated with 1 ML bismuth. For the H-termination, resulting in the formation of a Si(lll)-(lx 1)-H surface, 
a treatment in a 1% solution of HF acid is used, and after an additional cleaning step in deionized water the sample 
is transferred to the ultra-high vacuum within a time of 10 minutes. 

The distance-dependent four-point measurements are performed in a four-tip scanning tunneling microscope [S2) . 
In this system the individual positioning of the tips is realized under the control of an optical microscope. This has 
the advantage compared to an electron microscope that frequently observed influences of the electron beam on the 
surface properties are avoided |S31IS5| . In the linear arrangement of the four probes the voltage between the inner 
two tips is measured as a function of the current injected by the outer tips. For the square configuration, the voltage 
probes and current probes are located at the corners of opposite sides of the square, and by successively rotating the 
assignment of these probes different rotation angles in multiples of 90° can be realized for one fixed position of the 
tips. The four-point resistance is obtained from the slope of the measured I/V curves close to zero volt, whereby 
each resistance value is averaged over four I/V curves. 

As for the four-point measurements preferentially a non-equidistant probe spacing is used, in which three tips 
remain at a distance s and only the distance x between one outer current injecting tip and the adjacent voltage 
measuring tip is varied, Eq. 1 in the main article has to be modified in this case to |S6ItS8j 


with the 2D sensitivity X 2 D 


and the effective spacing s^j^ 


^20 = — -^ and = 


1 In 2 


i4p 


X2D T^<J2D 




(SI) 
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1.11 1 
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X 2s S -b X 


(52) 

(53) 


For contacting the sample surface, in a first step the four tips are approached until a tunnel contact with a current 
in the low nA range is established. Then, after retracting the tips by several nm, the feedback of the STM is switched 
off, and the tips are manually approached further separately until an increase in the current up to 1 /rA is observed. 
At this point, the tips are in contact with the sample, but they only touch the surface and have a penetration depth 
of only a few A as confirmed by separate experiments. For transport measurements, the sample is set to floating 
potential and a current is injected by an applied voltage ramp between the outer tips. Since the voltage at the inner 
tips is measured as function of the actual current, there is no influence of a potentially fluctuating contact resistance 
and resulting variation in current through the sample due to the switched-off feedback. 


ADDITIONAL EXPERIMENTAL RESULTS 

Additional four-point resistance measurements in the linear probe configuration are carried out on differently n- 
doped Si(lll)-(7x7) substrates with bulk resistivities ranging over two orders of magnitude. The results for the 
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FIG. SI. Measured four-point resistance of n-doped Si(lll)-(7x7) samples with bulk resistivities of (a) 20f2cm, (b) 2000 flcm 
and (c) 4100 flcm as function of the probe distances s and x for the equidistant conhguration as well as the non-equidistant 
configuration with s = 50 ^m. The decreasing four-point resistance in the equidistant region indicates a non-surface conductance 
contribution due to bulk and space charge layer. Applying a three-layer model for transport in the sample yields the solid 
red line located between the two limiting cases of a pure surface conductance (dashed blue line) and a pure bulk conductance 
(dashed green line). 
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non-equidistant probe spacing with a; < 50 /am and s = 50 /am, as well as the equidistant spacing with distances 
X = s > 50/am are shown in Fig. SI for three Si samples with resistivities of (a) 20ncm, (b) 2000 ffcm and (c) 


4100 ricm. From the best fits according to the used three-layer model (solid red line), the surface conductivities for 
the three different samples can be determined to tr/a) = (4.2 ± 0.6) • 10“® ct/;,) = (4.7± 0.6) • 10“® and 

CT(c) = (5.3 ± 0.8) • 10“® As the measured surface conductivity in the linear configuration is a combination of 

contributions due to step edges and the step-free surface, these values can differ slightly from each other for the three 
differently doped Si samples because of slightly different miscut angles of the substrates and resulting different step 
densities at the surfaces. Nevertheless, all obtained surface resistivities (including the 700 flcm sample discussed in 
the main text with a = (5.1 ± 0.7) • 10“® are very close to each other and still compatible within the error 

tolerances indicating that the step contributions are similiar for all samples. Comparing the measured data to the two 
limiting cases of a pure 2D conductivity with the above values and a pure 3D conductivity arising only from the bulk 
shown as dashed blue and green lines in Fig. |S1[ respectively, one can see that the observed resistance behavior is 
very close to the 2D case for all of the three differently doped Si(lll) samples, although the bulk resistance varies over 
several orders of magnitude. This indicates a separation of the surface layer from the bulk arising from the space charge 
region with low conductivity due to a depletion zone preventing an enhanced charge transport through the bulk, which 
especially becomes important for high bulk doping concentrations. For larger probe distances, the measured four-point 
resistance increasingly deviates from the 2D case and shows a pronounced non-surface contribution resulting from 
additional charge transport through space charge layer and bulk. This non-surface contribution decreases for lower 
bulk doping levels, as it is expected in comparison with the large increase in bulk resistance in the limiting case of a 
pure 3D conductivity. Within the three-layer model, the space charge region is approximated by one intermediate layer 
with thicknesses of (a) 0.9/rm, (b) 4.8/rm and (c) 5.3/rm, and with constant conductivities of (a) 1.5 • 10“'^ 

(b) 1.8-10“^ and (c) 3.2-10“"*^ respectively, for the three differently doped samples in Fig. SI A)-(c). 


THREE-LAYER MODEL FOR CONDUCTANCE 

The three-layer model assumes a layered sample structure shown in Fig. consisting of a thin surface layer, 
an intermediate layer and a semi-infinite bulk characterized by their respective conductivities cti, tT 2 and 173 , and 
positions of the interfaces Zi and Z 2 . At the surface a current I is injected by a cylindrical tip with radius r^. Due to 
calculation requirements, the surface layer cannot be two-dimensional, so that a finite thickness of one atomic layer 
(3 A) is assumed. As V • j = 0 for the current density j = ctE = —crV4) inside the sample (excluding the injection 
point), the electrical potential d* in this region can be determined by solving the Laplace equation 


A$ = 0 


(S4) 


in cylindrical coordinates. Taking account of the angle-independent polar symmetry for one tip, a solution for the 
potential in the individual layers is [S9] 
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with Jq denoting the Bessel function of the first kind. The corresponding boundary conditions are 


d 

cti^4>i(p,0) = -ji H{rt - p) 

(S8) 

<^i[p,zx) = <^2{p,zx) 

(S9) 

O-l^d*!)/?, Zl) = cr2^$2(p,2i) 

(SIO) 

^2{P,Z2) = (^■i{p,Z2) 

(Sll) 

0-2 J^d'2(p, ^ 2 ) = 0-3 J^$3(/?,Z2) 

(S12) 
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FIG. S2. The three-layer model consists of a thin surface layer, an intermediate layer and the semi-infinite bulk described by 
their respective conductivities cri, (T 2 and 0 - 3 , and their positions 21 and 22 . The current I is injected by a cylindrical tip of 
radius rt at the surface layer. 


resulting from the current injection (Eq. S 8 ), as well as from the continuous transitions of the potential (Eq. S9 and 
Eq. Sll) and the current density (Eq. SIO and Eq. S12) between the layers. In Eq. S 8 the expression H{rt — p) 
denotes the Heaviside step function. The injected current density can be described by ji = assuming a cylindrical 
tip with a tip radius of r* « 25 nm, which seems reasonable for an STM tip. Nevertheless, it turns out that also other 
values for the tip radius in the range of 5 nm to 100 nm do not influence the results of the calculations in a considerable 
manner. Based on Eqs. S 8 — S12 a matrix equation determining the coefficients a{k ),..., f{k) is derived 
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which can be solved by means of matrix inversion. 
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As the potential at the surface (2 = 0) can be expressed by 


(S13) 
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(S14) 


only the coefficients a{k) and b{k) are relevant for the calculation. From Eq. S13 an expression for the sum of these 
coefficients is obtained 
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Introducing cartesian coordinates with x = (x yf" and p = |x — Xo| = 1 /( 2 ; — Xq)^ + {y — yoY for a tip positioned at 
xq, the combined potential on the surface d>surf,i2 for a current source at position Xq^ and a current sink at position 
X 02 results in 


^’surf.l2(x) = $surf,l(|x - XqJ) - 4'surf.2(|x - Xq^ |) . (S 16 ) 

Finally, the four-point resistance measured on the surface is determined by the quotient of the potential difference 
between the positions X 03 and X 04 of the measuring tips, and the current I resulting in 


_ ^surf,12 (xOa ) 41 surf ,12 (XO4 ) 


(S 17 ) 
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FIG. S3. Description of the measured four-point resistance (n-doped Si(lll)-(7x 7) sample, bulk resistivity 700 Dcm) by the 
three-layer model for different input parameters (colored curves). For the nominal bulk conductivity and a space charge layer 
conductivity varying over three orders of magnitude from 1 ■ 10“® (magenta curve) to 1 • 10“® (green curve) 

the measured data are fitted with the surface conductivity as free paramter. All curves can describe the data in the region of 
small probe spacings (x < 50 /im) very well, while there are considerable deviations from the data in the region of larger probe 
distances (s > 50 ^m). Nevertheless, the obtained surface conductivity shows only a very minor spread and varies less than 
10% from (T 2 d = 4.8-10~® (magenta curve) to a 2 D = 5.4-10~® (green curve). This indicates that the measured 

four-point resistance in the small probe spacing region does not depend on the details of the space charge layer, and so the 
surface conductivity can be determined very precisely, even if the exact properties of the space charge layer are not known. 


For the linear probe configuration with equidistant spacing s between three tips and a non-equidistant spacing x 
between one outer current tip and the adjacent voltage measuring tip, Eq. |S18| simplifies to 


1 r°° 

i?'^^(s, x) = - [oi(fc) + b{k)] ■ [Jo(fcs) — Jo{k{s -I- x)) — Jo{2ks) + Jo(fcx)] dk . 

I Jo 


(S19) 


So, Eq. S19 describes the equidistant measurement range according to Fig. SI a)-(c), if x = s > 50 /im, and 


the non-equidistant region, if x < s = 50 /rm is assumed. The integral over the Bessel functions can be evaluated 
numerically and the result can be fitted to the measurement data with the free parameters CTi, CT 2 and Z 2 defining the 
properties of the surface layer and the intermediate layer, which is used as approximation for the space charge region 
in the sample. The value for cts is known from macroscopic measurements of the bulk resistivity and is in agreement 
with the nominal doping concentration. 

The surface conductance channel has the largest influence on the measured four-point resistance for small distances, 
i.e. X < 50 /im (region of non-equidistant probe spacing), while the properties of the intermediate layer (conductivity 
and thickness) have the largest influence on the four-point resistance in the region of larger distances (equidistant 
spacing region). This is illustrated in Fig. |S3[ in which the measured four-point resistance (n-doped Si(lll)-(7x7) 
sample, bulk resistivity 700 flcm) is compared to theoretical curves resulting from the three-layer model for different 
input parameters (colored curves). For the constant nominal bulk resisitivity of the sample and a varying conductivity 
for the intermediate layer (space charge layer) over three orders of magnitude from 1 • 10“^ (magenta curve) 

to 1 • 10“®f2“^m“^ (green curve), the measured data are fitted to determine the surface conductivity. If the space 
charge layer conductivity is enhanced above a value of 1 • 10“^ the measured data cannot be described any 

more by the model and a fit is not possible. So, the space charge layer should be described by a conductivity value 
in the chosen region. Although there is a very large spread in the space charge layer conductivity, the results show 
a very minor spread of the surface conductivity with a deviation of less than 10% from CT 2 d = 4.8 • 

(magenta curve) to a 2 D = 5.4 • 10“® (green curve). All fit curves describe the data very well in the region 

of small probe spacing (x < 50/tm), while there are considerable deviations from the data points for larger probe 
spacing (s > 50/tm). This confirms the influence of surface conductivity and space charge layer conductivity in 
different regions of the used probe distances. Thus the surface conductivity can be determined very precisely from 
the non-equidistant probe spacing region, even if the properties of the space charge layer are taken into account only 
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very approximately. On the other hand, the parameters of the space charge layer can be further approximated from 
the measurement data in the equidistant region. In total, this justifies the crude approximation of the space charge 
region as only one layer with constant conductivity, as the quantity to be determined, the surface conductivity cti, 
does not depend significantly on the values of the conductivity cr 2 and width Z 2 used to describe the intermediate 
layer. 

The three-layer model described above can easily be extended to a multi-layer model consisting of N separate layers. 
In this case, the boundary conditions in Eqs. |S8|- |SI2 have to be modified to include the transition between the layer 
n — I and layer n for n = 1,..., {N — 1) and the size of the matrix in Eq. S13 becomes {2N — 1) x (2iV — 1) requiring 
a numerical solution method, but Eq. |S18| for the four-point resistance remains the same. So, if a depth-dependent 
conductivity profile is known, the charge transport in a sample can be described very precisely. 

In contrast to the analytical model, one other way to simulate the contributions of the different charge transport 
channels to the four-point resistance is to invoke finite element calculations. We have not done this, since the distances 
involved range from a few nanometers (radius of the current injection) to 200 /im (maximum probe distance) and such 
a large range of length scales is difhcult to include in finite element calculations. Nevertheless, as the exact analytical 
solution of the potential problem in a layered sample has quite an elementary form, the three-layer model is easier to 
apply and provides a more accurate computation than a finite element simulation. 


* Corresponding author: b.voigtlaenderQfz-juelich.de 
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